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ABSTRACT 

RX J1713.7— 3946 is the most remarkable TeV 7-ray SNR which emits 7-rays in the highest energy 
range. We made a new combined analysis of CO and H I in the SNR and derived the total protons in 
the interstellar medium (ISM). We have found that the inclusion of the H I gas provides a significantly 
better spatial match between the TeV 7-rays and ISM protons than the H2 gas alone. In particular, 
the southeastern rim of the 7-ray shell has a counterpart only in the H I. The finding shows that the 
ISM proton distribution is consistent with the hadronic scenario that comic ray (CR) protons react 
with ISM protons to produce the 7-rays. This provides another step forward for the hadronic origin of 
the 7-rays by offering one of the necessary conditions missing in the previous hadronic interpretations. 
We argue that the highly inhomogcncous distribution of the ISM protons is crucial in the origin of 
the 7-rays. Most of the neutral gas was likely swept up by the stellar wind of an OB star prior to 
the SNe to form a low-density cavity and a swcpt-up dense wall. The cavity explains the low-density 
site where the diffusive shock acceleration of charged particles takes place with suppressed thermal 
X-rays, whereas the CR protons can reach the target protons in the wall to produce the 7-rays. The 
present finding allows us to estimate the total CR proton energy to be ^10"'® ergs, 0.1 % of the total 



energy of a SNe. 

Subject headings: cosmic rays - 

ISM: clouds ■ 

I. INTRODUCTION 



gamma rays: ISM ~ ISM: individual objects (RX J1713.7-3946) 
ISM: molecules - ISM: H I 



It is a long-standing question how the cosmic ray (CR) 
protons, the major constituent of cosmic rays, are ac- 
celerated in the interstellar space. Supernova remnants 
(SNRs) arc the most likely candidate for the acceleration 
because the high-speed shock waves offer a n ideal site 
for diffusive shock accelera tion (DSA) fe.g.. iBelll Il978l: 
iBlandford fc Ostrikeii [1971 . The principal site of CR 
proton acceleration is, however, not yet identified obscr- 
vationally in spite of a number of efforts to address this 
issue. 

RX J1713.7— 3946 is the brightest and most ener- 
getic TeV 7-ray SNR detected in the Galact ic plane 
survey with H.E.S.S. (jAharonian et al.l l2006a|) and is 
a primary candidate where the origin of the 7-rays 
may be established. Discovery of t he SNR was made 
in X -rays with ROSAT (Pfcffcr mann &: As chcnbacbl 
119961 ) and ASCA showed that the X-rays are non- 
thermal synchrotron emission with no thermal features 
(jKovama et al.l I1997D. TeV 7-rays were first detected 
by CANGAROO (jEnomoto et al.l [2001 and, subse- 
quently, H.E.S.S. resolved the shell-like TeV 7-ray 
distributi on with a ^ ? 1 degre e point spread function 
(JAharonian et al.l 12004 l2006bL |2007| ). The 7-rays 
are emitted via two mechanisms, either leptonic or 



hadronic, closely connected to the CR particles and it is 
important to understand which mechanism is working 
in the SNR. The leptonic process includes the inverse 
Compton effect of CR electrons which energize the low 
energy photons of the cosmic microwave background 
and additional soft photon fields (e.g., infrared). The 
hadronic process includes the neutral pion decay into 
7-rays following the reaction of CR protons with the low 
energy target protons in the interstellar medium (ISM). 
Considerable work has been devoted to explaining 
the 7-ray emission in the framework of leptonic and 
hadronic scenarios (lAharonian et al.ll2006bl: IPorter et al.l 



Ellison fc Vladimirov 2008; Tanaka et al 



2006t iKatz fc Waxman '2008' 'Bcrczhko fc VollJ 



2008; 



2008; 



Morlinoetall 120091: l Acero et al.l 120091: [Ellison et al 



201c ; Patnaude et al.in2010l: Zirakashvili fc Aharonian 
20TalAbdo et al.ll2011l:lFang et al.ll20M 

The molecular gas interacting with the SNR was dis- 
covered in Vlsr, the velocity with respect to the lo- 
cal standard of rest, around —7 km s^^ at 2.6 arcmin 
resolution based on NANTEN Galactic plane CO sur- 
vey and the distance of the SNR was determined to be 
1 kpc by using th e flat rotation curve of the Galaxy 
(jFukui et al.l I2003D . This determination offered a ro- 
bust verification for a small distance of 1 kpc and re- 
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vised an old valu e 6 kpc derived fr om lower resolution 
CO observations (jSlane et al.lll999l ). Studies of X-ray 
absorption suggested a similar distance 1 kpc under 
an assumption of u niform foreground gas distribution 
(jKovama et alJ[l997| ) , but the local bubble of H I located 
by chance toward the SNR makes the X-ray absorption 
uncertain in estimating the distance (jSlane et al.l 119991 : 
iMatsunaga et alj 120011) . A subsequent careful analy- 
sis of the X-ray absorption favors the smaller distance 
(jCassam-Chenai et al.ll2004| ). At 1 k pc the SNR has a 
radius of 9 pc and an age of 1600 yrs (jFukui et alJl2003l : 
iWang et al.lll997| ) and t he expanding shock front has a 
speed of 3000 km s~^ (iZirakashv ili fc Aharonian' 120071: 
lUchivama eFall I2003L |2007| ). ^loriguchi et al. (2005) 
showed further details of the NANTEN CO distribution 
and confirmed the identificati on of the inte r acting molec- 
ular gas at -20-0 km s~] - bv iFukui et~all (|2003f ). Most 
recently, iSano et aLl ()2010[) showed that the SNR harbors 
the star forming dense cloud core named peak C at sim- 
ilar T^jSR and argued that X-rays are bright around the 
core, reinforcing the association of the molecular gas. 

Importantly, the associated molecular gas with the 
SNR opened a possibility to identify target protons where 
the hadronic process is working. If the cosmic ray density 
is nearly uniform, we expect the 7-ray distribution mim- 
ics that of the interstellar target protons. Some nearby 
molecular clouds show good spatial correlation with rela- 
tively high resolution 7-ray images, clearly verifying that 
the hadronic process is working to produc e 7-rays in the 
cosmic ray sea (e.g., iBertsch et al.l [l993l and the refer- 
ences therein.). A detailed comparison between the ISM 
protons and the recent high-resolution 7-ray images of 
H.E.S.S. sources is useful to test the correlation, although 
such a test was not possible until recently in the preced- 
ing low r esolution 7-ray observatio ns at degree-scale res- 
olution. lAharonian et al.l (|2006bf ) compared the NAN- 
TEN CO distribution with the H.E.S.S. TeV 7-ray im- 
age and examined both leptonic and hadronic scenarios 
as the origin of the 7-rays. By adopting annular averag- 
ing of TeV 7-rays and C O in the shell (see Figure 17 of 
lAharonian et al.ll2006bD . these authors found that TeV 
7-rays are fairly well correlated with CO, whereas the 
correlation is not complete in the sense that the south- 
eastern rim of the TeV 7-ray shell has no counterpart in 
CO. The complete identification of target ISM protons 
thus remained unsettled in the hadronic scenario. 

The 7-rays and X-rays are significantly enhanced 
toward the clu mpy molecu la r gas at a pc scale as 
first shown by iFukui etlUI (|2003| ). A strong con- 
nection among the molecular gas, the 7-rays, the 
X-rays, and perha ps cosmic rays is therefore sug- 
gested (IFukui et al..-2003: Fuk^ i 2008; Sano et al.1 l2010l: 
IZirakashvili fc Aharonianll2010D . Most of the previous 
models of the 7- and X-rays cited earlier assume more or 
less uniform density distribution of the ISM in the SNR, 
whereas the observations of the ISM indicate that the 
actual distribution is highly inhomogeneous, with den- 
sity varying by a factor 100 or more around the SNR. 
In addition. Galactic-scale studies of 7-rays suggest that 
there is "dark gas" which is not detectable in CO or in 
H I b ut still contributes to the 7-rays ari d visual extinc- 
tion (|Grenier et al.ll2005l:lAde et al.l(201lD . Such gas may 
be either cold H I or II2 with no detectable CO. So, it 



is important to consider H I carefully in order to have 
a comprehensive understanding of the ISM protons. In 
the present paper, we shall use the term "dark H I" for 
observed H I with significantly lower brightness than the 
surroundings. So, "dark H I" does not necessarily mean 
"dark gas" above, while they might be linked. 

We here present a combined analysis of both the 
^^CO(J=1~0) and H I datasets in order to clarify the 
distribution of the ISM protons in RX J1713. 7-3946. 
The quantitative analysis is made mainly by using the 
^^CO(J=1-0) and H I data here and a comparative 
study of the ^^C0( J=l-0, 2-1) transitions is a subject 
of a forthcoming paper. A theoretical work of magneto- 
hydrodynamics which incorporates fully the inhomoge- 
neous ISM in RX J1713.7— 3946 is complimentary to the 
present observational paper and will be published sepa- 
rately (Inouc, Yamazaki, Inutsuka & Fukui 2011; here- 
after IYIF2011). The present paper is organized as fol- 
lows. Section [2] gives the description of the CO and H I 
datasets. Section |3| consists of five Sub-sections; Sub- 
section 3.1 gives a combined data set of CO and H I, 
Sub-sections 3.2-3.4 the distribution of total ISM pro- 
tons in the SNR with an emphasis on dark H I, the cold 
and dense atomic phase of ISM protons, and Sub-section 
3.5 comparison with the 7-rays. Section |4| consists of two 
Sub-sections; Sub-section 4.1 describes the initial distri- 
bution of the ISM before the stellar explosion, and Sub- 
section 4.2 the scheme of particle acceleration and its 
relationship with the 7-rays. The conclusions arc given 
in Section [5j 

2. DATASETS OF CO, H I AND TEV 7-RAYS 

The ^^CO(J=1-0) data at 2.6 mm wavelength were 
taken with NANTEN 4m telesco pe in 2003 April and ar e 
identical with those published bv lMoriguchi et aD ([2005') . 
The system temperature of the SIS receiver was ~250 K 
in the single side band including the atmosphere toward 
the zenith. The beam size of the telescope was 2f6 at 115 
GHz and we adopted a grid spacing of 2.0' in these obser- 
vations. We adopt hereafter the brightness temperature 
(K) as the spectral line intensity scale. The velocity res- 
olution and rms noise fluctuations are 0.65 km s^^ and 
0.3 K, respectively. 

The ^^C0(J=2-1) data at 1.3 mm wavelength were 
taken with NANTEN2 4m telescope in the period from 
August to No vember in 2 08 an d part of the dataset was 
published by iSano et aTl (j2010[ ) . The frontcnd was a 4 
K cooled Nb SIS mixer receiver and the singlc-side-band 
(SSB) system temperature was -^250 K, including the at- 
mosphere toward the zenith. The telescope had a beam 
size of 90" at 230 GHz. We used an acoustic optical spec- 
trometers (AOS) with 2048 channels having a bandwidth 
of 390 km s^^ and resolution per channel of 0.38 km s~^. 
Observations in ^■^C0( J=2-l) were carried out in the on- 
the-fly (OTF) mode, scanning with an integration time 
of 1.0 to 2.0 sec per point. The chopper wheel method 
was employed for the intensity calibration and the rms 
noise fluctuations were better than 0.66 and 0.51 K per 
channel with 1.0 and 2.0 sec integrations, respectively. 
An area of 2.25 square degrees in a region of 346.7 deg < 
/ < 348.2 deg and —1.3 deg < 6 < 0.4 deg was observed. 

The H I data at 21 cm wavelength are from 
the Southern Galactic Plane Survey (SGPS; 
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Fig. 1. — (a) The H.E.S.S. TeV 7 -ray distribution of RX J1713.7— 3946 in smoothed excess counts above the cosmic-ray baclcground (see 
Figure 2 of lAharonian et aLII2007lV Contours are plotted every 10 smoothed counts from 20 smoothed counts, (b) Averaged brightness 
temperature distribution of ^ CO(J=1-0) emission in a velocity range of Vlsr = —20 km s^^ to km s~^ is shown in color IjFukui et al.l 
120031 : iMoriguchi et al.ll2005l1 . White contours show the H.E.S.S. TeV 7-ray distribution and are plotted every 20 smoothed counts from 20 
smoothed counts, (c) Averaged brightness temperature distribution of H I emission obtained by ATCA and Parkes in a velocity range from 
^LSR = —8 km s~^ to —6 km s~^ IIMcClure-Griffiths et al.ll2005l ) is shown in color. White contours show the ^^CO(J=1-0) brightness 
temperature integrated in the same velocity range every 1.0 K km s~^ ('~3(t). 



iMcClure-Griffiths et all 120051 ) and combined from 
the Australia Telescope Compact Array and the Parkes 
Radio Telescope. The beam size of the dataset was 2.2 
arcmin and we adopted a grid spacing of 40" toward the 
RX J1713.7— 3946 in the current analysis. The velocity 
resolution and typical rms noise fluctuations were 0.82 
k m s~^ and 1.9 K , respec tively. 

IMoriguchi et alj (J2005D showed an analysis of the 
^^C0( J=l-0) distribution over 100 km s~^ with a coarse 
velocity window of 10 km s~^ in order to test association 
with the SNR. These authors showed that the velocity 
range Vlsr — —20 to km s^^ has convincing signs of 
association with the SNR. We adopt in the present work 
the velocity interval, —20-0 km s^^, for the associated 
ISM and present detailed ^^CO ( J=l-0 and 2-1) and H I 
data every 1 km s~^ (see Figure A in Appendix A). 

For the H.E.S.S. 7-ray data we us ed the combined 
H.E.S.S. image shown in Figure 2 of lAharonian et al.l 
(p007i) . Data of 2004 and 2005 are used for this 
smoothed, acceptance-corrected gamma-ray excess im- 
age. The TeV image utilizes minimum 3 H.E.S.S. tele- 
scopes in event reconstruction to obtain a Gaussian stan- 
dard deviation of 0?06 or FWHM of 0?14 (8.3'). 

3. COMBINED ANALYSIS OF THE CO AND H I 
DATA 

3.1. Distribution of CO and HI 

Figure 1(a) shows TeV 7-ray distribution toward RX 
J1713.7-3946 obtained by H.E.S.S. and Figure 1(b) 
shows a velocity averaged distribution of "'^^C0(J=1- 
0) overlayed on the TeV 7-ray distribution. The 
^^CO(J=1-0) intensity becomes larger in the north to 
the Galactic plane than in the south and the most promi- 
nent features above 0.7 K are located in the north- 
west. The general ^^CO(J=1-0) distribution is shell- 



like associated with the 7-ray shell, showing weaker or 
no CO emission in part of the south. There are two 
regions where -'^^CO(J=1-0) delineates particularly well 
the outer boundary of the shell in the southwest and east. 
In addition, we see some of the ^^C0( J=l-0) features are 
located within the shell including the prominent peak C 
at {I, h) = (347?07, -0?40). 

The ^^C0(J=2-1) distribution is qualitatively similar 
to the ^^C0( J=l-0) distribution. A typical ratio of the 
[J=:l-0]/[J=2-l] line intensities is ~0.6, consistent with 
what are derived in the other molecular clouds withou t 
heat source (e.g. lOhama et all 120101: iTorii et al.ll201lD . 
We tentatively choose from Figure A three major CO 
clouds, W, N and SW, and three minor ones, E, NE and 
SE, for the sake of discussion as schematically shown in 
Figure 2(a), where we use ^^C0(J=2-1) data by taking 
an advantage of higher angular resolution. 

Figure 1(c) shows an overlay of the H I distribution 
superposed on the ^^CO(J=1-0) intensity in a velocity 
range of —8.0 — 6.0 km s~^. The average H I brightness 
temperature ranges from 60 to 150 K and becomes higher 
toward the Galactic plane. The brightest H I of ^150 
K, the central cloud, is located toward the center of the 
SNR [(/, 6)=:(347?25, -0?38)] where little ^^GO{J=l~ 
0) is seen. We find dark H I clouds of around 60 K 
in the west (W cloud) and in the southeast (SE cloud). 
These dark H I clouds arc not due to absorption of the 
radio continuum radiation which is very weak toward the 
SNR (Lazendic et al. 200l). The dark H I W cloud weh 
corresponds to the ^^^00(7=1-0) distribution, showing 
sharp edges both toward the east and west. The dark 
H I SE cloud has almost no counterpart in CO. The 
relatively bright H I emission is seen in the north of the 
SNR (N cloud). The N cloud tends to be located toward 
^^C0( J=l-0) peaks, whereas the H I brightness shows a 
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Fig. 2. — (a) Schematic of the identified ^^CO(J=2-l) clouds is shown in colored contours. The image and white contours show the TeV 
7-ray distribution (Figure 1(a)). The integration velocity ranges are follows; —8 — 2 km s~^ (contour level: 2.3 K km s~^) for the SW 
cloud (CO), -17 — 5 km s"! (contour level: 4.9 K km s'^) for the W and N clouds (CO) and -6-0 km s"! (contour level: 2.7 K km s"!) 
for the NE, E and SE clouds (CO), (b) The locations of the identified H I clouds are shown in colored contours. The gray scale image and 
white contours show the TeV 7-ray distribution. The solid contours are for H I emission and the dashed contours are for dark H I. The 
integration velocity ranges are as follows; —8 — 2 km s~^ (contour level: 214 K km s~^) for the W cloud (H I), —12 — 4 km s~^ (contour 
level: 780 K km s"^) for the central and SE clouds (H I), and -14 — 11 km s'^ (contour level: 399 K km s'^) for the N cloud (H I). 



non-monotonic, more complicated behavior than in the 
W cloud. In the northeast, we find a rim of relatively 
lower H I brightness of ~100 K toward (/, &)=(347?5, 
— 0?25) that lies along the 7-ray shell. A schematic of 
the four main H I clouds is given in Figure 2(b). The 
good correspondence of the H I clouds with the CO and 
the 7-rays supports that the H I is physically associated 
with the SNR. 

In Figure 3 we show typical H I and CO profiles in the 
four main H I clouds. Figure 3 indicates that the H I 
emission is generally peaked at —10 km s~^ with small 
hints of saturation, confirming that the H I is associated 
with the SNR and is generally not optically thick. We 
find that narrow H I dips having depths of 20-30 K often 
correspond to ^^CO(J=1-0) emission features in the N 
and W clouds. The linewidths of the narrow H I dips 
are as small as a few km s~^. It is likely that these H I 
dips represent residual H I in cold CO gas seen as self- 
absorption. The broad H I dip in the SE cloud is also 
ascribed to self-absorption as argued into detail in Sub- 
section 3.3. We show H I expected profiles of the back- 
ground H I emission with a strai ght-line approximat ion 
as dashed areas in Figure 3 (c.g.. lSato fc Fukulll978f ). 

3.2. Molecular protons 

In order to covert the ^^CO(J=1-0) intensity into 
the total molecular column density, we use an X fac- 
tor, which is defined as X (cm^^/K km s^^) = A^(Il2) 
(cm~2)/W(^^C0) (K km s~^). In order to derive an 
X factor, the ^^CO(J=1-0) intensity is compared with 
the cloud dynamical mass (virial mass), or with the 7- 
rays produced via interaction of cosmic ray protons with 
molecular clouds. An X factor therefore accounts for the 
total hadronic mas s and is ob scrvationally uniform in the 
Galactic disk (e.g. iFukui fc Kawaniura i20m). We here 



adopt an X factor of 2.0x10^0 W{^^CO) (cm^^/K km 
s~^) derived from the 7-rays and -"^^C OfJ^^l-O) inten- 
sity in the Galaxy (jBertsch et al.lll993f ). We double the 
II2 column density to derive the ISM protons in molec- 
ular form as shown in Figure 7(a). Compared with the 
i2C0( J=l-0) line, the ^^CO (J=2-l) line is not a com- 
mon probe of the molecular mass. This is in part be- 
cause the "'^^C0(J=2-1) emission samples a smaller por- 
tion, having a higher excitation condition, of a molecular 
cloud than traced by the ^^C0( J=l-0) emission. We es- 
timate for instance that a typical fraction in area of the 
i^C0(J=2-l) emission to the ^'^CO{J=1~0) emission is 
about 70-80 % at the half-intensity level convolved to 
the same beam size in the present region from the CO 
data in Figure A. 

3.3. Atomic protons 
3.3.1. Optically thin case 

We use the 21cm H I transition to estimate the atomic 
proton column density. A usual assumption is that the 
H I emission is optically thin and the following relation- 
ship is used to calculate the H I column density; 



A^p(H I) = 1.823 X 10^^ I TL{V)dV (cm" 



(1) 



wher e Tl(V) is the observed H 1 brightness temperature 
(K) (jDickev fc Lockman|[l990l ). We note that this sim- 
ple assumption is usually valid and apply equation ([T]) 
to the regions where no H I dips are seen. It is certain 
that the narrow H I dips in the W and N clouds repre- 
sent self-absorption by cold residual H I in CO gas from 
their exact coincidence with CO in velocity. The most 
prominent dark H I cloud, the SE cloud, shows large 
linewidths, not so common as self-absorption. We shall 
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Fig. 3.— ^^CO{J=1-0) and H I profiles at the four H I clouds; the N cloud (l, b) = (347?50, -0?23), the SE cloud (l, b) = (347?57, 
-0?93), the central cloud (/, b) = (347? 27, -0?37), and the W cloud (l, b) = (347?07, -0?40). The positions are denoted in the H.E.S.S. 
TeV 7-ray distribution. The shaded area shows expected profiles behind the self-absorption. 
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Fig. 4.— (a)The H.E.S.S. TeV 7-ray distribution toward the SE cloud ifXh 



aronian et al. '2007l ) . Red contours show averaged H I brightness 
temperature distribution in a velocity range from —15 km s~^ to —5 km s~^ (McClur c-Griffiths et aOl2005l ). (b) The H I and ^^CO(J=1-0) 
spectra at {I, b) =(347?55, — 0?92). The shaded area shows an expected H I profile. 



examine if the SE cloud represents self-absorption in the 
followings. 

3.3.2. The dark H I SE cloud 

We first show the integrated intensity image of the SE 
cloud in Figure 4(a). The H I contours are every 3.9 a 
noise level and shows significant details not apparent in 
Figure 1(c), where a coarser color code is used. We find 
the H I brightness variation is generally well correlated 
with the shell of 7-rays in gray scale in the sense that H I 
brightness decreases toward the enhanced 7-rays. This 
trend lends a support for physical connection of the SE 
cloud with the 7-ray shell and may be interpreted as due 



to decrease in spin temperature with density increase in 
the self-absorbing H I gas (see Sub-section 3.3.3). No 
^^CO(J=1~-0) emission is seen toward the SE cloud, ex- 
pect for a possible small counterpart at (Z, b) = (347?64, 
-0?72) and Vlsr = -6-0 km 8"^ (Figure 2(a) and Figure 
A), suggesting that density of the SE cloud is lower than 
the CO clouds. 

Figure 4(b) shows a typical H I profile in the SE cloud 
having a deep and broad dip. The large velocity span 
of 20 km s~^ is not so common as a self-absorption 
feature; in nearby dark clouds H I self-absorption is 
generally narrow with a few km s~^ in linewidth (e.g. 
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Fig, 5, — (a)j4v distribution l|Dobashi et aLfeOOSI ) is shown in color 
density of the total ISM protons estimated from both CO and H I in a velocity range from Vlsr 
self-absorption is talsen into account. Contours are the same as in Figure 1(a). 



Contours are the same as in Figure 1(a). (b) Distribution of column 

-20 to 10 km s-^ Here the H I 



iKrco fc Goldsinithll2010[) . wliereas H I self-absorption as 
broad as 10 km s~M s seen in giant molecular clouds (e.g. 
iSato fc Fuk"mlll978[) . The SE cloud delineates the 7-ray 
shell (Figure 4(a)) and is possibly compressed gas by the 
wind of a high-mass star, the SN projenitor. We have 
investigated the velocity distribution of the SE cloud as 
given in Appendix A. We find that the SE cloud shows a 
strong velocity gradient which matches the blue-shifted 
part of an expanding swcpt-up shell. Such a shell is a 
natural outcome of the stellar-wind compression by the 
SN projenitor, supporting that the broad H I dip is as- 
cribed to the acceleration of H I gas by the wind. A H I 
stellar-wind shell in Pegasus driven by an early B star 
indeed show s a linewidth as large as 15 km s^^ (Sub- 
section 4.1, lYamamoto et al.j I2OO60 . similar to that of 
the SE cloud. The difference from the narrow H I dips 
in the W and N clouds may be due to density; the SE 
cloud has lower density and is subject to stronger accel- 
eration than the CO clouds with narrow H I dips (see 
for further discussion Sub-section 4.1), whereas the CO 
clouds having higher density arc less accelerated by the 
wind, making a systematic velocity gradient less clear in 
CO than in H I (sec Figure B3). 

Figure 5(a) shows the distr ibution of the extin ction 
Ay toward RX J1713. 7-3946 (jDobashi et al.l[2005l ). and 
indicates that the SE cloud, as well as the rest of the 
shell, is traced by the enhanced optical extinction. This 
lends another support for the self-absorption interpreta- 
tion of the SE cloud. Figure 5(b) shows the total (molec- 
ular and atomic) ISM proton column density both in 
the SNR (derived later in Sub-section 3.4) and in the 
foreground within 1 kpc, which is supposed to corre- 
spond mainly to the optical extinction. The total pro- 
ton column density A^p of ^ 10^^ cm~^ in Figure 5(b) 
corresponds to extinction of ^4 magnitude if we adopt 
the relationship A^n(cm~ ^) = 2.5x10^^ ■ Ay (magnitude) 
(j Jenkins fc Savaga[ll974l ). The extinction toward the SE 
cloud is 2-3 magnitude in Figure 5(a) and is consistent 



with the H I self-absorption by considering the contami- 
nation by the foreground stars which tends to reduce Ay 
toward the Galactic plane. 

In summary, we find it a reasonable interpretation that 
the SE cloud represents H I self-absorption associated 
with the SNR shell. 



3.3.3. Analysis of the HI self-absorption dips 

We shall briefly review some basic properties of H I 
gas in or der to understand th e behavior of H I bright- 
ness (e.g. iSato fc Fukiulll978D . The spin temperature, 
Tg, of H I is ^100 K or higher in warm neutral medium 
at particle density less than 10 cm~'^. Tg decreases with 
density from 100 K down to 10 K in a density range o f 
100-1000 cm-3 (e.g.. Figure 2 in lGoldsmith et al.l[2007[ ). 
The temperature decrease is mainly due to higher shield- 
ing of stellar radiation and increased line cooling. 

It is well established that H I is converted into II2 
on dust surfaces with increasing of the gas column den- 
sity and UV shielding and that II2 is dissociated by cos - 
mic rays and UV photons (e.g.. lAllen fc Robinsonl[l977| ). 
The equilibrium H I abundance is determined by the bal- 
ance between formation and destruction of H2 and the 
residual density of H I is about 10~^ that of H ^ in typi- 
cal interstellar mole cular clouds (jAllen fc Robin son 1973 
ISato fc Fukulll978| ). We also note that the II2 abun- 
dance should be time dependent since the formation of 
II2 is a slow process in the order of 10 IVlyrs fo r density 
around 100 cm'^ (^e.g.. lAllen fc Robinson[|1977[) . 

Based on the H I-H2 transition, we interpret the dark 
H I in Figure 2(b) as representing the H I with lower 
Tg- The CO W cloud shows a good spatial coincidence 
with the dark H I W cloud as is consistent with the 
interpretation. The other prominent dark H I region, the 
SE cloud, shows no CO and we suggest that its density 
is lower and its Tg is higher than in the CO W cloud. H I 
brightness Tl{V) is expressed as follows (e.g., Sato and 
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Fig. 6. — (a) Distribution of peak optical depth of the H I self-absorption, (b) Distribution of atomic proton column density, Afp(H I), 
estimated for the H I emission and sel f-absorption. The veloc ity range in the both figures is from —20 km s~^ to km s~^. Contours show 
the H.E.S.S. TeV 7-rays distribution l lAharonian et aT]|2007l ) and are plotted at 20 smoothed counts. We assume spin temperatures Ts of 
40 K and 10 K, inside and outside the dotted box toward the SE cloud, respectively 



Fuliui 1978); 

+ [rf G(T/) + rBG]e-^(^) - (T<EG + T^G) (2) 

where Tl(F), T^, t{V), TJG(^) and T^'^iV) are the 
observed H I brightness temperature, the spin temper- 
ature, the optical deptli of cold H I in the cloud, and 
the foreground and background H I brightness temper- 
ature, respectively, at velocity V. T^ and T^'^ arc 
the continuum brightness temperature at 21-cm wave- 
length in the foreground and background of the cloud, 
respectively. The r adio continuum emis sion is weak in 
RX J1713.7-3946 (iLazendic et al.l[200l . and T^° and 
Tq'^ are nearly zero as compared with Tl{V). 

We are then able to estimate the H I column density 
of dark H I clouds. Figure 4(b) shows the H I self- 
absorption dip with the background H I emission inter- 
polated by a straight line connecting the two shoulders 
at and 20 km s~^. This gives a conservative estimate 
because the actual background H I shape perhaps has a 
more intense peak at —10 km s~^ as seen in the northern 
area of the SNR. The spin temperature Tg of the dark H I 
gas is an unknown parameter. We estimate Tg to be less 
than ~55 K from the lowest H I brightness at the bottom 
of the dip in Figure 4(b) and higher than ~2 K, where 
the t emperature of the CO clouds is ^10 K (|Sano et al.l 
|201C1[ ). We estimate the absorbing dark H I column den- 
sity to be A'^p(H I) = 1.0 xlO^^ cm"^ (optical depth = 
0.8), 1.8x1021 cm-2 (optical depth = 1.1) and 3.1x10^1 
cm~2 (optical depth =1.5) for assumed three cases Tg 
= 30, 40 and 50 K, respectively, for the half-power line 
width Aw=10 km s"^, where the H I optical depth f is 
estimated by equation (2) and A^p(H I) by the following 
relationship; 

A^p(H I) (cm^2) ^ i_g23 x lO^^T^ (K) Au(km s-i)f(3) 



We shall here adopt Tg = 40 K and a corresponding dark 
H I optical depth of 1.1. A higher Tg gives a higher 
optical depth and vice versa. The relatively large optical 
depth around 1 is consistent with the fairly fiat H I dip 
in Figure 4(b), which suggests weak saturation. We also 
tested the effects of elevating the background H I by 15 
K and found a small change of 5 x 10^" cm"^. The error 
is mainly introduced by the straight-line approximation 
and uncertainty in Tg of -^10 K. We infer the dark H I 
column density is accurate within a systematic error of 
~lxl02i cm-2. 

The average H I density in the SE cloud is roughly esti- 
mated to be 150 cm^^ by dividing 1.8 x 10^^ cm"^ by ^4 
pc, the line of sight length of the thick ISM shell, follow- 
ing the three-dimensional model described in Sub-section 
3.5. This density is significantly lower than the critical 
density for collisional excitation of the CO ( J=l-0) tran- 
sition, ^1000 cm~'^, consistent with no CO emission from 
the SE cloud and with low spin temperature around 40 
K. 

We also extended such an analysis to the regions with 
narrow H I dips associated with CO emission, where we 
adopt Ts = 10 K, the kinetic temperature of the CO 
gas. The small dips in these regions indicate that the 
H I optical depth is generally as low as ^0.1 reflecting a 
small fraction of the residual H I in CO gas. We show 
the distributions of the peak optical depth of the H I 
self-absorption in Figure 6(a), and the derived total H I 
column density distribution, the sum of the H I in emis- 
sion and self-absorption in Figure 6(b), where the SE 
cloud is significant. We shall hereafter refer to the dark 
H I of Tg = 40 K as "cool H I" and that of Tg = 10 K as 
"cold HI". 

3.4. Total ISM protons 

The number of the total ISM protons in the SNR is 
given by summing up the three components in a velocity 
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Fig. 7. — (a) Distributions of column density of ISM protons A''p estimated from -'^^CO(J=1-0) Afp(H2), (b) H I emission with correction 
for the H I self-absorption A'^p(H I) and (c) sum of Np(R2) and Afp(H I). All the datascts used here are smoothed to a HPBW of TeV 
7-ray distribution with a Gaussian function, (d) TeV 7-ray distribution. Contours are plotted every 50 smoothed counts from 20 smoothed 
counts. 



range from —20 to km s^^; H2 derived from ^^C0( J=l- 
0), dark H I (dips) and warm H I (emissions). The re- 
sults are shown as spatial distributions in Figure 7. Fig- 
ure 7(a)-(d) show iVp(H2), iVp(H I), Np{}l2+il I) and 
TeV 7-rays, respectively. We see the total ISM protons 
.^p(H2+H I) shows a shell-like shape similar to the TeV 
7-rays which significantly improves the correlation with 
the 7-rays as compared with the case of molecular gas 
only. We therefore conclude that the contribution of H I 
is critical as well as H2 in counting the ISM protons. We 
find that in the south the total ISM proton is dominated 
by the atomic gas, whereas in the north the molecular 
and atomic protons are both important. A more quan- 
titative comparison will be given in Sub-section 3.5.2. 
Similar diagrams of the total ISM protons to Figure 7 



are presented for the optically-thin case for reference Fig- 
ure CI in Appendix C, where the shell-like distribution 
toward the SE cloud is missing. 

3.5. The j-rays and the ISM protons 
3.5.1. Gamma-ray distribution 

The TeV 7-ray distribution obtained by H.E.S.S. is 
a nearly circular-symmetric shell with some ellipticity 
elongated in the north-south direction. In order to gain 
an insight into the distribution of the 7-ray emissivity 
we undertake a simple analysis of the 7-ray distribution. 
We fi rst adopt an ellip t ical ann ular ring in the analysis, 
while lAharonian et al.l ()2006bD made a similar analysis 
by using a circular annular ring in correlating 7-rays and 
NANTEN CO intensity (see their Figure 17). 



22.0 



a 




347.5 347.0 

Galactic Longitude [degree] 



700 



600 



I 500 



.| 400 

o 
o 
E 
^300 



200 

100 


-1 



-•- TeV y-ray 

-B- Np(m 

^ iVp(Hi) 

T^ Afp(H2+Hl) 



M 



3*r 






-^ 



-A-^-^^k--^ . ^v-^^_ 



^-: 



-0- 



-ffi- 



.■^: 



-B 



(b): 



>* 



i-B- 



14 

12' 

10 

8 

6 

4 

2 





(U 

Q 
c 

E 

o 
O 

c 
o 

o 



80 -90 90 

Azimuth Angle [deg.] 



inside ring 
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Fig. 9. — Radial distribution of TeV 7-rays radiation. Small dots 
show the distributions of all the H.E.S.S. data points and large 
filled circles with error bars show averaged values at each radius. 
We assume a 3D spherical shell with a Gaussian-like intensity dis- 
tribution along its radius to approximate the TeV 7-ray distribu- 
tion (see text). The green line shows the estimated 3D Gaussian 
distribution and the red line shows its projected distribution. The 
peak radius and the full width at half maximum of the green line 
are estimated to be 0.46 deg (~ 8.0 pc) and 0.24 deg (~ 4.2 pc), 
respectively. 



We estimated the radius of the 7-ray shell as defined at 
a half-intensity level of the peak 7-ray smoothed count 
every 15 degrees for an assumed center. We averaged 
the radii in angle and minimized the sum of the squares 
of the deviation from the average. This process gives a 
central position to be (/, b) = (347?34, -0?52). For this 
central position, we plotted the radius every 15 degrees 
and found that a sinusoidal distribution is a reasonable 
approximation as expected. Fitting this plot by a si- 
nusoidal curve, wc find the shell is approximated by an 



elliptical shape with an aspect ratio of 1.1 whose major 
axis is almost in the north-south direction. This elliptical 
shape is adopted in Figure 8(a). 

Figure 9 shows the radial scatter of 7-ray smoothed 
counts and an averaged value shown by a step function 
in radius r every 0.05 degrees. Here we also adopted 
the elliptical shape and normalized the radius to that 
of the major axis with the elliptical modification. After 
several trials of different functional forms, we found a 
Gaussian radial distribution of the 7-ray emissivity per 
volume reproduces well the projected radial distribution 
in Figure 9. In the fitting we have two free parameters 
of the Gaussian shape, the peak radius tq and the sigma 
a expressed as follows; 



F(r) = ylxe-('^-''«)'/2"' 



(4) 

where A is a normalization coefficient. By requiring that 
the error in the fitting becomes minimum in the pro- 
jected distribution shown by the step function, we found 
ro = 0.46 degrees and cr = 0.10 degrees give the best fit 
as shown in Figure 9. This distribution shows that the 
observed shell is consistent with a shell of a half-intensity 
thickness '^0.24 degrees with nearly zero emissivity to- 
ward the center. This analysis indicates that the 7-rays 
are mainly emitted in a thick shell of 8.0 pc radius and 
4.2 pc width at the half-intensity level with nearly zero 
emission from the inn er part. A similar thick-s hell model 
was also obtained bv lAharonian et all ()2006b[) . Numeri- 
cal modeling of the 7-ray emission has been undertaken 
by several authors and indicates that the 7-ray emis- 
sion has a rather steep gradient beyond the peak of the 
shell in cither of the l eptonic or hadronic scenario (e.; 
iJun fc Norinaiil 119961 : IZirakashvili fc AharonianI I201C 
The fitting to the H.E.S.S. data above shows that the 
gradient in the 7-ray distribution is not so steep toward 
the outside, which may be due to smearing in space by 
averaging. We shall not try a further elaborated analysis 
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Fig. 10. — Radial distributions of averaged values of TeV 7-rays 
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show column densities estimated from ^^CO(J=l— 0) and H I, re- 
spectively, and 7Vp(H2+H I) shows the total ISM proton column 
density. 



here due to the Umiting angular resolution of H.E.S.S. 
which is 0.14 deg (FWHM). 

Figure 9 shows that the projected radial distribution 
of ISM protons follows a fairly similar distribution to the 
7-rays inside the SNR. This is consistent with that the 
ISM distribution is also shell-like with an inner cavity 
as is consistent with the stellar wind shell discussed in 
Sub-section 4.1; if the ISM has no cavity in the inner 
part, the projected distribution of the ISM should in- 
crease toward the center. We shall assume hereafter that 
the ISM distribution is also approximated by the same 
Gaussian shape as the 7-rays with a radius of 8.0 pc with 
a thickness of 4.2 pc at the half-intensity level. 

3.5.2. Comparison between the -y-rays and the ISM protons 

In the hadronic scenario, the target distribution should 
be correlated with the 7-ray distribution for a uniform 
CR distribution. This correlation should be seen inside 
the shell of the SN shock which has a sharp gradient be- 
yond its outer radius. We expect that the ISM protons 
are distributed beyond the outermost edge of the shell 
where CR protons cannot reach by diffusion. Beyond 
the SNR shock, the 7-ray emission profile may be influ- 
enced by components from the diffuse cosmic-ray back- 
ground and by the energy dependent transport of escap- 
ing cosmic -rays from RX J17 1 3.7— 3946 into the clum py 
ISM (e.g.. IGabiei ei~all[2009l : rCasanova et ani2010al lbh. 
We are able to avoid possible effects of such cut-off by 
taking the radius of the correlation analysis well within 
the SNR shell where the CR protons do not decrease in 
energy density. 

The ISM proton distribution is shown in Figure 8(a) 
with the two annular elliptical rings along the shell, 
where the size of the outer ring was chosen to meet the 
requirement above. Figure 8(b) shows a comparison be- 
tween the ISM protons and 7-rays in the position angle 
shown in Figure 8(a), where the vertical scale is adjusted 
so that the correspondence with the TeV 7-rays becomes 
optimum. Here, the error in the TeV 7-ray emission from 
the publicly available H.E.S.S. image is approximately 
(smoothed counts)°'^. In Figure 8(b) the uncertainty in 
the dark H I in the SE cloud, 1x10^^ cm~^, is in the or- 



der of 10-20 % of the total. The total ISM proton density 
shows a good agreement with the TeV 7-ray angular dis- 
tribution and also the central part in the inner ring. We 
recall that CO alone showed marked deficiency toward 
the S E cloud as com pared with the 7-rays (see Figure 
17 of lAharonian et al. 2006bl ). The present analysis in- 
dicates the deficiency is recovered by including H I and 
has shown that the total gas of both atomic and molec- 
ular components have a good correlation with the TeV 
7-rays in the annular ring. The total mass of the ISM 
protons responsible for the 7-rays is 2.0 xlO"^ Mq over 
the whole SNR (radius 0.65 deg); the mass of molecu- 
lar protons is 0.9x10^ Mq and that of atomic protons 
is 1.1x10''' Mq, where we assume the ISM protons inter- 
acting with the CR protons is proportional to the TeV 
7-rays (Sub-section 3.5.1., Figure 10). 

There arc two points in Figure 8(b), for which addi- 
tional remarks may be appropriate. One is the point at 
an azimuth angle of 115 degrees which may be estimated 
too low due to lack of correction for the self-absorption 
because of the large velocity shift in the expanding shell 
(see Figure Bl). Another is the point at an azimuth an- 
gle of 1 65 degrees where the strong CO emission (peak 
A after iFukui et al.l |2003[ ) increases the proton column 
density, although the increased protons may not be in- 
teracting with the CR protons beyond the SNR shock, 
leading to less 7-rays. 

An independent test is made by the radial distribution 
of the ISM protons given in Figure 10, where an average 
taken over the same binning as the 7-rays in Figure 9 
is shown by a step function and the total ISM protons 
and 7-rays are superposed with the same proportional 
factor as adopted in Figure 8. Here, the error in the TcV 
7-ray emission is approximately (oversampling-corrected 
total smoothed count)°'^ normalized to 1 (arcmin)^. We 
see the Np(R2+ii I) and 7-rays show a good agreement 
inside the shell and the 7-rays sharply decrease outside 
the shell. This offers another presentation of the good 
correlation between the 7-rays and the ISM protons. 

We argue that the apparent anti-correlation between 
the H I brightness at the bottom of the dips and the 7- 
rays in the SE cloud (Figure 4) is consistent with that 
the H I dips are due to the cool and dense H I gas. The 
anti-correlation is interpreted that the spin temperature 
Tg of H I decreases with density (Sub-section 3.3) and 
that the 7-rays increase with the ISM proton density lo- 
cally in the SE cloud, demonstrating detailed correspon- 
dence between the 7-rays and the ISM protons which is 
mainly atomic. The small and narrow H I dips in the W 
and N clouds have the H I column density less than 10^° 
cm^^, significantly lower than the typical molecular col- 
umn density by two orders of magnitude. So, in most of 
the regions except for the SE cloud the H I column den- 
sity is dominated by emission but not by self-absorption. 
For the sake of reference, we show a set of similar dia- 
grams of ISM proton distributions for the optically-thin 
case in Figures C2 and C3 in Appendix C, corresponding 
to Figures 8 and 10, respectively. 

Before concluding this Sub-section, we cautiously note 
that the cool/cold H I could not be estimated accurately, 
if the cool/cold H I is optically thick, if the cool/cold 
H I lies behind optically thick foreground H I in the line 
of sight, or if the background H I profile has a differ- 
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TABLE 1 
A Comparison between R,X J1713. 7-3946 and Pegasus Loop 



RX J1713.7-3946* Pegasus Loopt 



Distance (kpc) 1 

Diameter (pc) 17.4 

Total mass of the ISM (M©) ~20000* 

Thickness of the ISM shell (pc) ~4.2 

Peak brightness of H I (K) ~170 

Lincwidth of H I (km s"!) ~20 

Expansion velocity the gaseous shell (km s~^) ~10 

Spectral type of the projenitor Bl V / BO V* 



0.1 

25 

~1500 

~5 

~40 

~16 

-7-9 

B2 IV 



Note. — * ll''ukui et aEI ||2003| '). IMoriguchi et al.l ||2005|') . t lYamamoto et al.l (12003 '). ** 
ICassam-Chenai et al.l II2004I ) , ^ present wor k. The Pegasus loop ma y consist of two shells and the the 
mass should be regarded as an upper limit IIYamamoto et aI1l2006l) . 



ent shape from its neighbors. Such effects, while posing 
intrinsic hmits for probing cool/cold H I, are relatively 
unimportant for nearby objects at a distance of 1 kpc or 
less where foreground H I is not important. The dark H I 
W and SE clouds are probably good examples where the 
cool/cold H I is well traced by the low H I brightness, 
whereas the N cloud with higher H I brightness may be 
partially affected by the foreground H I in the line of 
sight. 

4. DISCUSSION 
4.1. The evacuated cavity by the stellar wind 

It is likely that the CO shell in Figure 1(b) was formed 
over a timescale of Myrs by the stellar wind of the pro- 
jenitor, an OB star that exploded as a supernova (SN) 
1600 yr ago. The total velocity span of the CO shell, 
^20 km s^^, is much smaller than the SN shock speed 
and indicates that it takes Myr to form the shell of the 
ISM as roughly estimated by dividing the radius 9 pc by 
10 km s~^. Molecular gas expanding at 10 km s~^ can 
move only 0.01 pc in 1000 yrs. Therefore, the current CO 
distribution has little been affected by the supernova ex- 
plosion (SNe) and holds the initial condition before the 
shock interaction. 

While a stellar-wind shell with a known central star 
is not often observed elsewhere, one such example is the 
Pegasus loop found in ^^C0( J=1~0), H I and dust emis- 
sion at {I, b) = ( 109°, —45°) cent ered on a run-way star 
HD886 (B2 IV) ()Yamamoto et a l. 2006). The Pegasus 
loop is located at ~100 pc in a relatively uncontami- 
nated environment outside the Galactic plane. No SNe 
occurred yet in this shell. A comparison between RX 
J1713.7— 3946 and the Pegasus loop is given in Table 1. 
In Pegasus the swept-up shell of the ISM has a width of 
~5 pc for a radius of ~18 pc and a total mass of ^1500 
Mq . The shell is mostly atomic and con sists of 78 smaller 
J^CO(J=1-0) clumps (see Figure 10 in lYamamoto et aTl 
|2006[ ). The clumped CO is a natural outcome of ther- 
mal/gravitational instability and seems common in such 
a shell. The shell is expanding at a total velocity span 
of 15 km s~^. The H I density inside the shell is ^1 
cm~^ in the north, where the stellar wind evacuated the 
ISM over 1 Myrs. The Pegasus loop is located in a some- 
what lower-density environment than RX J1713.7— 3946 
and offers an insight into the initial condition of the ISM 
prior to the SNe in RX J1713. 7-3946. 



Inoue, Yamazaki & Inutsuka (|2009D and IYIF2011 car- 
ried out numerical simulations of the hydro-dynamical 
interaction between the shock wave and the highly in- 
homogeneous neutral gas to model the interaction in 
RX J1713. 7-3946. The SN in RX J1713.7-3946 ex- 
ploded in the cavity with average density less than 1 
cm~5 fe.g.lZirakashv ili fc A haronian|[20Tol : IMorlino et al.l 
[20091 IBerezhko fcVolk 2Q0|) and the dense shell with 
CO clumps remaining more or less as they were prior 
to the SNe. The SN shock front moves almost freely 
at > 3000 km s~^ in the cavity in the early phase of 
^1000 yrs and begins to interact with the dense and 
thick clumpy ISM wall swept-up by the stellar wind only 
in the last few 100 yrs. The 7-ray shell is not strongly 
deformed, while we see some deviations of a pc scale 
from a perfect circular shell, suggesting effects of recent 
dynamical interaction. 

The interaction between molecular clumps and the 
shock is observed as the X-ray enhancement around 
dense molec ular clumps at a spatial resolution higher 
than 0.5 pc. iSano et all (|2010l ) showed that the molecu- 
lar clump peak C is rim-brightened in X-rays, suggesting 
that it is a dense clump overtaken by the shock, and 
peak A (jFukui et al.ll2003l) is also X-ray brightened only 
toward its inner edge, indicating the shock interaction 
at the inner boundary of peak A. IYIF2011 showed that 
the initial magnetic field i? of 1 fj,G is amplified to 0.1 
to 1 niG near dense clumps by the enhanced turbulence 
driven by the shock. The stronger magnetic field ex- 
plains the X-ray enhancement as due to the enhanced 
synchrotron emission that is proportional to i?^, or, due 
to increased acceleration. IYIF2011 also showed that 
the shock speed Vg is significantly reduced locally with 
density n (cm~^) such that Vs ~3000 km s~^ /y^n/no, 
where uq =1 cm~^. This dependence of Vg on density 
can explain the absence of thermal X-rays in the SNR 
because the molecular gas is too dense to be affected by 
the shock to emit thermal X-rays (IYIF2011). A uni- 
form lower-density case with si gnificant thernial X- rays 
by shock heating is presented bv lEllison et al.l ()2010[ ) but 
such a model is not applicable to the highly inhomoge- 
neous ISM of RX J171 3.7-3946 (IYIF201 1, see also dis- 
cussion in Section 4 of lEllison et al.|[2"010( ). The picture 
above is also consistent with that peak C, having den- 
sity greater than 10^ cm"^, has survived without erosion 
(ISano et al.l[20Tol) . 
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4.2. The j-ray emission mechanism 

TcV 7-rays arc emitted via two mechanisms, either 
leptonic or hadronic processes. The leptonic process 
explains 7-rays via the inverse Compton effect between 
CR electrons and low energy photons. In the hadronic 
scenario 7-rays are emitted by the decay of neutral 
pions which are produced in the high energy reactions 
between CR protons and ISM protons. Diffusive shock 
acceleration (DSA) is t he most w idely accepted scheme 
of pa rt icle acceleration (lBelllll978 : Blandford & Ostrikeij 
[19781 iJones fc EIIIsmJ Il991t iMalkov fc Drurv. .20011 
The previous works on RX J1713.7— 3946 show that 
the observed spectral energy distribution of 7-rays and 
X-rays is explained by either of the leptonic and/or 
hadronic mech anisms if DSA works t o accelerate 
the particles (TAharoiiian et al.l l2006bt iPorter et al.l 



2OOI iKatz fc WaxmanI I2OOS 
Ellison fc Vladimirovl 12005 ' 



IBerezhko fc Vo!^ 



Tanaka et al.l 



2008 : 



2008' 



Morlino et al.l 120091: lAcero et al.l 120091 lEllison et al 



2010: Patnaude et al. 2010; Zirakashv ili fc Aharonian 
~010 : Abdo et al. 2011 : Fang et al. 2011) 



In the hadronic scenario, where the neutral pion decay 
determines the 7-rays via proton-proton reactions, the 
average density of the target protons is constrained by 
the total energy of CR protons; the average target den- 
sity greater than 0.1 cni^"^ is required to produce CR pro- 
tons having the total energy of 10^^ erg, for the maximum 
energy of a SNe, while higher target density is required 
for less CR proton energy. In the leptonic scenario, where 
the inverse Compton process produces 7-rays, the critical 
parameter is the magnetic field which constrains the syn- 
chrotron loss timescale of CR electrons; a magnetic field 
of or der of 10 /iG is usually required (e.g. iTanaka et al.l 
[2008h . 

We here argue that the highly inhomogeneous dis- 
tribution of the ISM, the cavity and the dense and 
clumpy wall opens a possibility to accommodate the 
low-density site for DSA and the high-density target si- 
multaneously as discussed into detail by IYIF2011. A 
similar argument on the hadronic interaction between 
CR protons with the ambient dense c louds has been 
presented by I Zirakashvili fc Aharoniaiil (|20100 . In this 
picture, first, the cosmic rays are accelerated via DSA 
in the low density cavity, and second, the CR protons 
reach and react with the target protons in the dense 
wall to produce 7-rays. The main energy range of the 
CR prot ons required for hadronic TeV 7-rays is 10- 
800 TeV (jZirakashvih fc Aharonia^l2010[ ). The penetra- 
tion depth, Zpd, of cosmic rays is expressed as follows 
(IYIF2011); 



'pd 



O.It] 



1/2 



10 TeV J 



B 



100 ^G 



-1/2 



/ lag 



1/2 



V103 yr, 

(pc) (5) 



where E, B and iago are the particle energy, the magnetic 
field and the age of the SNR. The parameter 77 is the 
so-called " gyro- factor" and has some ambiguity. In the 
SNR, it i s reasonable to conside r 77 '^1 at least around 
the cloud lUchivama et al.l ()2007f ). Thus, the penetration 
depth of the protons in the above energy range is 0.3-2.8 
pc for magnetic field of 10 /iC and 0.1-0.9 pc for 100 /iG 
in a typical timescale of ^10"^ yr. The penetration depth 



of the CR electrons is determined by taking tagc equal to 
the synchrotron loss timescale (e.g. iTanaka et al.l 120081 ) 
in equation (6) and becomes energy-independent for the 
X-ray emitting electrons of 1-40 TeV as follows; 



/ = 0.026 7] 



1/2 



B 



100 AfG 



-3/2 



(pc) 



(6) 



We estimate I to be from 0.8 pc for 10 /iC to 0.026 pc 
for 100 /iC if 77 = 1. CR protons can therefore reach 
and penetrate into the dense gas within pc-scale of the 
acceleration site to produce TeV 7-rays, while the CR 
electrons stay relatively closer to the acceleration site, 
in particular, near the dense gas having strong magnetic 
field. This offers an explanation on the hadronic 7-ray 
production and the correlation between the 7-rays and 
target protons in Figures 4, 8 and 10 is a natural outcome 
i n the scenario (IYIF 2011). 

iGabici et al.l ()2007f ) discussed the importance of the 
energy-dependent inte raction between CR protons and 
molecular clouds and I Zirakashvili fc Aharonian! (|2010f ) 
discussed that the 7-ray spectrum may not distinguish 
the leptonic and hadronic scenarios in case of RX 
J1713.7— 3946 due to such energy dependence. Recently, 
Fermi-LAT observations showed that the GeV spectrum 
of RX J1713.7— 3946 is hard, simi lar to what is ex - 
pected in the leptonic scenario, and lAbdo et al.l (|201lD 
discussed that the hard spectrum may favor to the lep- 
tonic scenario. IYIF2011, however, argued that the hard 
Fermi-LAT GeV spectrum is explained well also by the 
hadronic scenario as due to the energy-dependent pene- 
tration of CR protons into the dense clouds and that the 
leptonic scenario is not unique to explain the spectrum. 
IYIF2011 confirmed that the 7-ray spectrum becomes 
similar both for the leptonic and hadronic scenarios, 
not usable to distinguish t he tw o scenarios, as noted by 
I Zirakashvili fc Aharonian! (|2010( ) and concluded that the 
hadronic origin is testable only by comparing 7-rays with 
the ISM target distribution. The present results have 
demonstrated that the ISM proton distribution show in- 
deed a good spatial correspondence with the 7-rays by 
taking into account the contribution of the H I and match 
with the prediction by IZirakashvili fc Aharonian! ()2010l ) 
and IYIF2011. 

The total energy of CR protons is estimated by the re- 
lationship between the total target protons and the ob- 
served 7-rays (2-400 TeV) afte r extrapolating the proto n 
spectrum to 1 GeV as follows (jAharonian et al.]l2006b( ) : 



Wt, 



1 - 3 X 10 



50 



Ikpc 



1 cm 



(erg) (7) 



where the distance to the source d ~1 kpc and the density 
of the target protons is n. The average density of ISM 
protons is calculated to be '--'130 cm~^ for the total mass 
of the ISM protons 2.0 x lO" Mq over the whole SNR (ra- 
dius 0.65 degrees) as modeled in Figure 10 and the total 
CR proton energy to be ^^0.8-2.3 x 10** ergs by using 
equation ^. This corresponds to ^ 0.1 % of the total 
energy release of a SNe and may appear low. The other 
SNRs like W44 and W28 of a few to 10 times 1000 yrs old 
have the t otal CR proton energy in the order to lO''^- 
10^° ergs (lAbdoet all [2010! : IGiulianieFaniMoh . We 
may speculate that the CR protons become accumulated 
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in a few times 10000 yrs to reach more than 10 % of the 
SNc energy. This issue is to be further tes ted by exam- 
ining cosmic ray es caping from SNRs (e.g. iGabici et al.l 
I200a lCasanovaet~a l. 2010a.bl). 

To summarize the discussion, we have shown that a 
combined analysis of CO and H I provides a reasonable 
candidate for the target ISM protons and thereby lends 
a new support for the hadronic scenario. We should note 
that the present analysis offers one of the necessary con- 
ditions for the hadronic scenario for uniform CR pro- 
ton distribution, but it is not a full verification of the 
hadronic scenario and does not rule out leptonic compo- 
nents. We need to acquire additional observations before 
fully establishing the hadronic scenario, including better 
determination of the magnetic field and higher angular 
resolution images of 7-rays at least comparable to that 
of the ISM. Chercnkov Telescope Array will provide such 
images in future. We discussed that the observed highly 
inhomogeneous distribution of the ISM plays an essen- 
tial role in the 7-ray production; DSA works in highly 
evacuated cavity and the accelerated CR protons travel 
over a pc to interact with the surrounding dense ISM 
protons. It is important to develop a similar analysis of 
both H I and CO in the other similar objects like RX 
J0852.0-4622 (Vela Jr.), RCW86 and HESS J1731-347. 
Such works are in progress based on the NANTEN2 ob- 
servations and high-resolution H I interfcrometry. 

5. CONCLUSIONS 

We summarize the main conclusions as follows; 

1. A new analysis of CO and H I has revealed that 
the TeV 7-ray SNR RX J1713. 7-3946 is associated 
with a significant amount of H I gas without H2 
derived from CO. This H I gas is relatively dense 
and cold and detectable mainly as H I emission. We 
have also identified regions where H I is observed 
as dark H I in self- absorption dips and derived the 
total ISM proton column density over the SNR. 
The H I plus II2, the total ISM protons, provides 
one of the necessary conditions, target protons, in 
the hadronic origin of the 7-rays. Such target ISM 
protons have not been identified in the previous 
study that took into account only II2, although the 
present finding alone does not exclude the leptonic 
origin. 

2. For an annular pattern around the TeV 7-ray shell, 
we compared the total ISM proton distribution 
with the TeV 7-ray distribution and found that 
they show reasonably good correspondence, vary- 
ing by similar factors. The inclusion of the atomic 
protons observed as the H I self-absorption dips is 
essential particularly in the southeast of the 7-ray 
shell. The interpretation of H I self-absorption dips 
is also supported by the enhanced optical extinc- 
tion toward the southeast rim. 



3. The cavity surrounding the SNR was created by the 
stellar wind of the SN projenitor. The inside of the 
cavity is of low density with <1 cm""^ while the cav- 
ity wall consists of the dense and clumpy atomic or 
molecular target protons of > 100-1000 cm^''. The 
diffusive shock acceleration in the highly inhomo- 
geneous ISM offers a reasonable mechanism of par- 
ticle acceleration in the low-density cavity and the 
dense wall acts as the target for 7-ray production 
by the CR protons. Hydro-dynamical numerical 
simulations of the interaction have shown detailed 
physical processes involved (IYIF2011). 

4. By considering the other pieces of the observational 
and theoretical works accumulated thus far, the 
present results make the hadronic interpretation 
much more comfortable in RX J1713.7— 3946. The 
current energy of the total CR protons is estimated 
to be ~10''^ ergs, 0.1 % of the total energy of SNe, 
if we assume the 7-rays arc all produced by the 
hadronic process. 
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APPENDIX 

APPENDIX A 
VELOCITY CHANNEL DISTRIBUTIONS IN RX J1713. 7-3946 

We show velocity channel distributions of "'^^CO(J=1-0, 2-1) and H I every 1 km s~^ from —20 km s^^ to km s~^ 
superposed on the TeV 7-ray distribution in Figure A. 



14 



FUKUI ET AL. 



0.0 



-0.5 



-1.0 



O -^ 

ATCA&Parkes NANTEN NANTEN2 H.E.S.S. 4pc 

•20.0- -19.0 km s' -19.0- -18.0 kms' -18. 0- -17.0 km s' 



Hi Brightness Temp. pC] CO Brightness Temp. [K] 



55 110 165 0.0 2.5 5.0 7.5 
17.0- -16.0 km s' -16.0- -15.0 km s' 



.0 O 



c- 



o 



o 



d- 



o 



=0 



^ 



& 



C3 



o 



^ 



■9 



o 

% 










347.5 347.0 
Galactic Longitude [degree] 



347.5 



347.0 



347.5 



347.0 



347.5 



347.0 



347.5 



347.0 



Fig. a. — Velocity channel distributions of ^■^CO(J=1-0, 2-1) emission and H I brightness temperature overlaycd on the TeV 7-ray 
distribution. First-row panels (top): H I image and '^■^CO(J=l— 0) contours. Second-row panels: H I image superposed on the TeV 
7-ray contours. Third-row panels: ^■^CO(J=1—0) image superposed on the TeV 7-ray contours. Fourth-row panels: ^■^CO(J=2-l) image 
superposed on the TeV 7-ray contours. Each panel shows CO and H I distributions every 1 km s~^ in a velocity range from —20 to km 
s~^. The lowest contour levels of CO and TcV 7-rays arc 0.73 K (~ Scr) and 20 smoothed counts and contour intervals of CO and TeV 
7-rays are 0.73 K (~ 3(t) and 10 smoothed counts, respectively. 



APPENDIX B 
EXPANDING MOTION OF THE DARK H I SE CLOUD 

Figure Bl left shows schematically an expanding spherical shell of radius i?o = 9 pc and uniform expansion velocity 
Vb = 10 km s~^ and Figure Bl right a position- velocity diagram of the shell, where the ellipsoidal nature of the shell is 
not taken into account for simplicity. Figure B2 shows three representative velocity-channel distributions of the dark 
H I SE cloud for a velocity range from —20 to —10 km s~^ and shows that the SE cloud is extended to the north. 
The extension shifts toward the northwest with velocity decrease from —10 to —20 km s~^ as is consistent with the 
iso- velocity contours expected from the shell model in Figure Bl. Figure B3 shows another presentation of kinematical 
details of the SE cloud in position-velocity diagrams. We choose a line AB passing through the center of the SNR 
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and the SE cloud, and another line CD passing through the SE cloud in the north-south (Figure B3(a)). We show a 
position- velocity distribution of H I along the line AB (Figure B3(b)) and H I profiles along the two lines AB and CD 
(Figure B3(c)). We find the SE cloud is extended to the northwest with a large velocity gradient of 10 km s^^ per 
0.5 degrees, or ~1.2 km s^^ pc^^. The H I profiles in Figure B3(c) shows that the dips are deep and clear at b less 
than —0.5 degrees but becomes shallower above b = —0.5 degrees. The shallower dips make it nontrivial to quantify 
the dips at b higher than —0.5 degrees; we note that, even when the dips are not clearly seen, the H I probably suffers 
from self- absorption to some extent as suggested by the weaker H I brightness at —12 km s~^ toward b = — 0?52 than 
toward b = — 0?35 (line AB). We note that the strong velocity gradient in Figures B2 and B3 is consistent with the 
blue-shifted part of an expanding shell. The strong velocity gradient is interpreted in terms of the expanding shell as 
depicted by a white circle in the position-velocity diagram (Figure B3(b)). The blue shift by 10 km s~^ toward the 
center of the SNR indicates this part of the shell is in the foreground. This is consistent with that the dips are due to 
self-absorption against the background H I emission. We also infer that the swept-up shell is highly non-uniform since 
the broad H I dips are seen only in a quarter of the shell. 
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APPENDIX C 
ANALYSIS OF THE H I EMISSION; THE OPTICALLY THIN CASE 

The present analysis has shown that the H I is self-absorbed in part of the SNR as indicated by the H I dips and the 
H I column density is estimated by taking into account the self-absorption (Figure 7). In order to see the effects of 
the self-absorption quantitatively, we here show for comparison the ISM proton distribution in the optically thin case, 
which does not take into account the self- absorption. Figure CI, equivalent to the self- absorption case in Figure 7, 
includes the H I column density distribution for the optically-thin assumption smoothed to the HESS resolution ((b) 
and (c)), where the SE cloud is not seen. Figures Cl(a) and (d) are the same with those in Figure 7. Figure C2 is 
equivalent to Figure 8. Figure C2(a) is the total ISM proton column density for the optically thin H I at NANTEN 
resolution overlayed on the TeV 7-ray distribution. Figure C2(b) is the corresponding azimuthal distribution of ISM 
protons and TeV 7-rays, where the ISM protons is deficient in azimuthal angle from —90 to degrees as compared to 
Figure 8(b). Figure C3 is equivalent to Figure 10, and shows the radial distribution of ISM protons for the optically thin 
H I without correction for the H I self-absorption. In the smoothed radial distribution, the effect of the self-absorption 
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is not so obvious. 
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Fig. B1.- 
a radius of 



— Schematic image of a uniformly expanding shell and its velocity distribution in the position-velocity plane. Here we assume 
the shell Rq and an expansion velocity Vcxp of 9 pc and 10 km s~^, respectively. 
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Fig. B2. — (Left) Velocity channel distributions of H I integrated intensity toward the SE cloud superposed on the TeV 7-ray contours. 
TeV 7-ray contours are plotted every 10 smoothed counts from 20 smoothed counts. The faded area in the upper panel is a component 
unrelated to the SNR. (Right) Model velocity distributions of an expanding shell shown in Figure Bl. Iso-velocity lines are shown here, 
and blue areas show the corresponding velocity range shown in the left panels. 
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Fig. B3. — (a) Aver aged brightness t emperature distribution of H I in a velocity range from —20 km s^^ to km s~^. Contours show the 
H.E.S.S. TeV 7-rays I IAharonian et aT1 ^2007') and are plotted every 10 smoothed counts from 10 smoothed counts. The line AB is inclined 
by 60 degrees to the Galactic plane and the line CD passes the center of the SNR. (b) Position-velocity distribution of H I along the line 
AB in Figure B3{a). The velocity resolution is smoothed to 1 km s^^ and the integration interval is 200 arcsoc. The white circle shows a 
schematic image of an expanding spherical shell (Figure Bl). (c) H I spectra along the lines AB and CD in Figure B3(a). Expected profiles 
of H I self-absorption are shown by straight lines in the spectra with significant FI I dips. 
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Fig. CI. — (a) Distributions of column density of the ISM protons A^p estimated from ^^CO{J=l— 0) Np(li2), (b) H I emission without 
correction for the self-absorption Afp(H I) and (c) sum of A'^p(H2) and Afp(H I). Here we assume for reference that the H I emission is 
optically thin and the H I self-absorption is not taken into account. All datasets used here are smoothed to a HPBW of the TeV 7-ray 
distribution with a Gaussian function, (d) TeV 7-ray distribution. Contours are plotted every 50 smoothed counts from 20 smoothed 
counts. 
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Fig. C2. — (a) Distribution of column density of ISM protons Afp(H2+H I) in a velocity range from —20 Ism s^^ to km s^"', where 
the H I is assumed to be optically thin and without self-absorption. Contours and two elliptical rings are the same as in Figure 8(a). (b) 
Azimuthal distributions of A'^p(H2), Np{}i I), A'^p(H2-|-H I) and TeV 7-ray smoothed counts per beam in the two elliptical rings in Figure 
C2(a). The same plots inside of the inner ring arc shown on the right side in Figure C2(b). 
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Fig. C3. — Radial distributions of averaged values of TeV 7-rays, A^p{H2), A'^p{H I) and A'^p(H2-|-H I), where the H I is assumed to 
be optically thin as in Figure CI. Afp(H2) and Np(H I) show column densities estimated from ^^CO(J=1-0) and H I, respectively, and 
Wp(H2+H I) shows the total ISM column density, the sum of A''p(H2) and Afp{H I). 



